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Abstract. 

Optical vortices generically arise when optical beams are combined. Recently, we 
reported how several laser beams containing optical vortices could be combined to form 
optical vortex loops, links and knots embedded in a light beam (Leach et al 2004 Nature 
432 165). Here, we describe in detail the experiments in which vortex loops form 
these structures. The experimental construction follows a theoretical model originally 
proposed by Berry and Dennis, and the beams are synthesised using a programmable 
spatial light modulator and imaged using a CCD camera. 

PACS numbers: 03.65.Vf, 02.10.Kn, 42.25.Hz, 42.40.Jv 

Optical vortices (phase singularities, wave dislocations) are ubiquitous in nature, 
occurring whenever three or more light waves interfere. They are places in the 
interference field where the intensity is zero, and the phase is undefined, and usually 
occur at points in 2- dimensional fields, and along lines in 3 dimensions (see, for 
example, El Ej)- Within the laboratory, vortices are commonly produced using 
computer generated holograms or diffractive optics. Holograms of the design originally 
demonstrated by jHIH] are now widely used to create light beams with vortices embedded, 
and the orbital angular momentum [HI IB] with which they are associated. In 2001, 
Berry and Dennis showed theoretically 01101 that specific superpositions of beams could 
be generated in which vortex loops could be linked together and even knotted. Here, 
we describe in detail our experimental generation ^T] of such laser beams. 

An optical vortex embedded on the axis of a laser beam is usually characterised 
by a function, R\ l \ exp(iZ0) describing the local structure of the optical field around the 
vortex point in a transverse cross section of the beam. The exp(iZ</>) factor describes 
the phase singularity nature of the vortex, B)" shows that this must be accompanied 
by a zero of intensity. The integer / determines the number of phase cycles made in a 
right-handed circuit of the vortex and is called the strength of the singularity. Although 
I can take any integer value, positive or negative, vortices with \l\ > 1 are unstable 
and, in the presence of astigmatism, degenerate to form multiple vortices of I = ±1 
upon propagation. Many observations have been reported on the seeming interaction 
between two or more vortex points as light propagates from one plane to the next (see, for 
instance, ^2] and references therein). In this view, two vortices with oppositely signed 
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/ may spontaneously appear together in one plane, propagate to a subsequent plane, 
and then recombine (birth/creation and death/annihilation of vortices). By identifying 
propagation of the beam with temporal evolution, this language is possibly misleading 
since the optical beam is static in time in the laboratory reference frame. 

A different physical description is to recognise that the pair of vortices forms a 
looped thread within the beam (see, for example, |13j). An arbitrary light field, such 
as a space-filling speckle pattern, contains many optical vortex lines in a complicated 
arrangement jHj . As parameters vary, the vortex lines evolve in a complicated way, and 
may reconnect, and loops may appear or disappear [T| ITUl [T5[ ITS]. 

The original analysis of knotted and linked optical vortex loops [Qj used a 
combination of Bessel beams JTj. From an experimental standpoint this construction 
has two problems: firstly, that ideal Bessel beams have infinite aperture (their transverse 
profile, specified by a J/ Bessel function, is infinitely extended); secondly, that the range 
of intensity in the proposed beams was too large to observe the faint light around the 
knotted and linked dark threads. The construction was later extended (12] to give the 
same topological structures within paraxial beams, allowing the possibility of knotted 
dark structures within Laguerre-Gaussian beams ;6j, with fixed frequency specified by 
the wavenumber k. In cylindrical coordinates R, 0, z, these have the normalised form 



, , p\ jgexp(#) ( -R 
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tt(|/| + p)\ {w 2 + iz/kf^ 1 * \2{w 2 + iz/k) 
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where L denotes the generalized Laguerre polynomials. In addition to an axial optical 
vortex of strength I, there are two other parameters that may be varied - the radial 
number p, which is an integer determining the order of the Laguerre polynomial 
(effectively, the number of nodal rings about the beam axis), and the waist width w, 
which is a continuous variable setting the finite transverse width of the beam. These 
additional degrees of freedom allow us to create the same topological structures, but 
with a reduced range of intensity, enabling the nodal structures to be unambiguously 
seen. 

We synthesize the link and knot by superposing three Laguerre-Gaussian beams 
with the same index I, giving the field ^unpert, then perturbed by another Gaussian 
beam ?/> pert with index I = 0. Therefore, the knot and link fields V'knot.iink are equal to 
the sum V'unpert+V'pert- The amplitudes, values of p and waist width w of the components 
of the unperturbed superposition V'unpert are chosen so that at a specific radius Rq in 
the waist plane, both the field and its radial derivative are zero, i.e. 

V'unpert(-Ro) = 0, <9i<:?/\i n p er t|R = 0. (2) 

The resulting destructive interference forms a configuration of three circular vortex 
loops in the transverse plane, centred on the z-axis. One is in the waist plane 
(z = 0) with radius Rq, and the other two, which share the same smaller radius, are 
symmetrically arranged on either side of this plane. The fact that the radial derivative 
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vanishes on the central vortex ring (z = 0, R = Rq) means that it has zero strength - 
it can be thought of as two superposed rings with opposite strengths. The numerical 
value of Ro and particular choices of p and w in the superposed beams are chosen such 
that the intensity is relatively bright in the regions separating the three vortex rings 
from each other and from the /-fold singularity on the axis. 

To this superposition is added a simple Gaussian perturbing beam V'pert with near 
planar phasefronts (i.e. wide waist). As the amplitude of the perturbing beam increases, 
the configuration of three vortex rings deforms through a sequence of reconnections 
(explained in detail in [TU]) to give an (1,2) torus knot [18J. Therefore, the total 
superposition of beams we used is 



The superposition of beams with amplitudes ai,2,3 form VVpert, satisfying the conditions 
(J2J) above. For convenience, ai = 1, and wq is a fixed transverse lengthscale of the beam, 
chosen in our experiment to be 1mm. The link, depicted later in figures |U and El was 
formed with I = 2, and the relative amplitudes: a 2 = —0.35, a 3 = 0.36, a pert = 0.35 
(which gave Rq = 2wq). The knot of figures El and was formed with I = 3, and relative 
amplitudes a 2 = — 0.26, a 3 = 0.29, a pert = 0.25 (which gave R = 2.2w ). If I = 1 
the 'torus knot' is a single loop, wrapped twice around the vortex axis, which will not 
be described further. The beam combinations we have chosen, satisfying the general 
conditions of (TU], are not unique (two equations are solved, but three parameters (a 2 , a 3 
and Rq) are varied; also, the values of p and w can be changed). They are convenient in 
that they are straightforward to implement and the resulting intensity range makes it 
possible to image the singularities. Adjustment of the superposition numerically reveals 
that that the link occurs for values of the relative perturbation amplitude a per t between 
0.3 and 0.42. Both the knot and link are stable between similar finite ranges of their 
various parameters. 

As discussed above, vortex carrying beams can be easily created within the 
laboratory using diffractive optical components. For example, if a diffraction grating 
is modified such that the lines near its centre combine to form a fork dislocation, then 
the first-order diffracted beam contains an optical vortex. The strength I of the vortex 
is equal to the order of the fork dislocation. In addition, p circular discontinuity or 
discontinuities can be added to these designs to give the radial index of the resulting 
beam [T9*] . Although a hologram of this kind provides precise control over the phase 
of the diffracted beam, the intensity distribution depends also on the profile of the 
illumination beam. It transpires that for a simple single-forked hologram design, 
illumination by a fundamental Gaussian mode gives an I = 1 vortex beam in the far field 
which has over 80% of its energy in the I = l,p — Laguerre- Gaussian mode. Although 
satisfactory for many studies and applications, this level of mode purity is insufficient 
to form the topological structure we desire. 



^link,knot( a; ; V i z ) ~ ^unpert + ^pert 



ai^ /0 (x, y, z; 2w ) + a 2 ipio{x, y, z; w ) 




(3) 



Vortex Knots in Light 



4 



The design of phase-only holograms can be modified to control not only the phase 
structure of the diffracted beams but their intensity also I2U 122]- in essence the 
intensity at any point in a plane can be attenuated by adjusting the efficiency of the 
blazing in that region of the hologram. 

Our required superposition of Laguerre- Gaussian beams could be generated with 
multiple holograms, but it is more convenient to produce a single diffractive component 
that produces the combination of beams directly. We adopt a simple yet highly 
effective approach to phase-only hologram design. Once the required indices, amplitude 
weightings and relative phases of the ingredient modes are determined theoretically (as 
described above), the phase distribution of the link and knot beam in the waist plane, 
&(x, y)beam = & r g ^iink,knot ( x , V ■, 0) , is calculated. To this phase distribution is added the 
phase distribution of a blazed diffraction grating A) gra ti ng (where A is the period 
of the grating), such that the first order diffracted energy is angularly separated from 
the other orders, which are subsequently blocked using a spatial filter. The desired 
intensity I(x, y) of the superposition is also calculated and normalised so the maximum 
intensity is unity. This intensity distribution is applied as a multiplicative mask to the 
phase distribution of the hologram, acting as a selective beam attenuator imposing the 
necessary intensity distribution on the first order diffracted beam. 

The resulting phase distribution of the hologram y)hoio> expressed within a 
to 2ir range is given by 

Hx,y) hol0 = ((($(x,y) bcam + A) grating ) mod2 ^. - tt)) 

x sine 2 ((1 - I(x,y))n) + n. (4) 

The sine distortion of the intensity terms accounts for the mapping of phase depth to 
diffraction efficiency. An example of this process is shown in figure ^ 

To implement this hologram design we use a programmable spatial light modulator 
(HoloEye LC-R3000) under computer control to produce the desired beam. Such devices 
are based on liquid crystal arrays, which introduce spatially dependent phase retardation 
across the beam profile, allowing the creation of any specific combination of beams. 
Compared to a static hologram, the spatial light modulator allows us to control the 
beam parameters dynamically and, most importantly, the amplitude of the perturbation 
in real time. 

The optical beams were produced using a spatial filtered He-Ne laser, expanded to 
overfill the 20mm aperture of the spatial light modulator, and configured to generate 
the appropriate superposition of Laguerre-Gaussian beams, as shown in figure El 

The hologram was calculated at 800 by 800 pixels and displayed on the SLM at 
1200 by 1200 pixels. The period of the diffraction grating was 0.4 lines/mm such that, 
after passing through a lens, the first order diffracted energy was angularly separated 
from the other orders by approximately 2mm at the Fourier plane of the hologram. The 
unwanted diffraction orders were subsequently blocked by centering the first diffracted 
order on an adjustable iris which acted as a spatial filter. The size of the spatial filter was 
such that diffraction effects were negligible. A second lens was used to create an image 
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of the modulator. A CCD array positioned in this second plane and translated along the 
beam axis measured the intensity structure of the beam on each side of the waist. Given 
that the fine details of the vortex threads lie in regions of near darkness and in order 
to measure the vortex positions within the beam cross-section, it was necessary to over- 
saturate the CCD. The resulting image was mainly white, but contains points of darkness 
corresponding to the vortices as they intersected the plane of the CCD. A comparison 
of experimental and theoretical plots demonstrating the increasing saturation of the 
detector in the waist plane of the link field are shown in figure El Also shown here is the 
measured phase pattern of the beam, showing that the dark vortex threads are indeed 
singularities of the optical phase (this phase distribution was measured according to the 
methods described in 

These transverse vortex coordinates were measured in neighbouring planes, enabling 
the 3-dimensional spatial configuration to be deduced. The amplitude of the 
perturbation was chosen initially at the theoretically optimized value, then adjusted 
to give the best defined link or knot. In figures E] and El images taken at various 
transverse planes through the vortex link and knot beams are shown alongside their 
simulated counterparts. The measured positions of the vortices are extremely close to 
those determined theoretically for the modelled beams. 

Figures El and [7| represent the experimental spatial vortex link and knot 
constructions reconstructed from the measurements. Represented in this way, it is 
clear that the vortex topologies are as claimed, although the vortex loops making them 
are not particularly smooth. This is a consequence of the choice of superposed beams, 
here chosen so that the range of intensity in the vicinity of the vortex loops is small 
enough to pick out the vortices. In other constructions, such as that in [9 j using Bessel 
beams, the knot and link curves are smoother, but sit in regions of near-darkness. It 
is possible that smooth vortex knots may exist in constructions similar to ours with a 
reduced range of intensity, but we were unable to find these. 

We have described the creation and observation of static knots and links of dark 
vortex threads in light waves, verifying an earlier prediction of Berry and Dennis that this 
was possible within optical wave fields. Knotted vortices have been studied theoretically 
since Lord Kelvin's vortex atom hypothesis j23] in a range of different physical situations, 
such as hydrodynamics (23120], field theory [21] and nonlinear excitable media I29j. 
However, in these contexts, nonlinearities mean that vortex knots are probably unstable. 
Our structures are superpositions of free space optical modes with the same frequency, 
so are temporally stable. It is not known theoretically whether arbitrary topological 
configurations of vortices can be embedded in free space optical fields; only a few special 
cases are known (the knots and links studied here 131101, which have to be threaded, and 
certain braids J3U], which require counterpropagating beams). In addition to being of 
fundamental interest, these experiments illustrate topological light shaping that might 
be used as a topological waveguide for quantum mechanical matter waves, such as Bose- 
Einstein condensates jHU E2] • 
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Figure 1. Illustration of the algorithm used to generate the holograms used in the 
link experiment (the knot hologram was similar). The desired phase of the beam 
at the beam waist is added mod27r to a carrier. This is then multiplied by the 
desired intensity of the beam (corrected to account for the mapping of phase depth to 
diffraction efficiency) to produce the hologram which is used in the experiment. The 
0-2-7T range is represented in greyscale, as is actually programmed into the hologram. 
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Figure 2. Experimental configuration used to produce the vortex beams. The 
expanded and collimated beam from a He-Ne laser is reflected by the spatial light 
modulator displaying the required hologram. The beam is passed through a lens to 
spatially separate the diffraction orders. The first diffracted order is selected, and the 
other orders removed, by passing the light through an adjustable iris. A second lens 
is used to recollimate the beam and allow transverse cross-sections of the beam to be 
imaged. A CCD camera is placed at different planes to record the intensity of the 
beam. 
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Figure 3. Experimental and theoretical representation of the transverse structure 
of the waist of the link field. A sequence of plots of intensity, with increasing degree 
of saturation, shows how the only structure measured are the vortex lines piercing 
the observation plane. The phase (conventionally represented using hue), shows the 
cyclic phase change around the vortices, establishing the singular phase nature of the 
intensity nodes. 
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Figure 4. Theoretical and experimental transverse cross-sections of the vortex link. 
The distances refer to the distance from the beam waist. In the model, the positions 
of the vortices are indicated by coloured spots, red (I = +1) and blue (I — —1). These 
vortices are observed as black dots in the saturated model and the experimental results. 
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Figure 5. Theoretical and experimental transverse cross-sections of the vortex knot. 
The distances refer to the distance from the beam waist. In the model, the positions 
of the vortices are indicated by coloured spots, red (I = +1) and blue (I — —1). These 
vortices are observed as black dots in the saturated model and the experimental results. 
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Figure 6. 3-dimcnsional representation of the experimental vortex link, reconstructed 
from the measured data. Different projections of the link are shown, and the direction 
of propagation in each case is indicated by the arrow. The link is formed by the red and 
green loops. The threading vortices, which do not take part in the link, are represented 
by thinner blue and grey lines. The scale for the direction of propagation has been 
reduced by approximately 1000. 
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Figure 7. 3-dimensional representation of the experimental vortex trefoil knot, 
reconstructed from the measured data. Different projection of the knot are shown, 
and the direction of propagation in each case is indicated by the arrow. The knot is 
represented by the red loop. The threading vortices, which do not take part in the 
knot, are represented by thinner blue and grey lines. The scale for the direction of 
propagation has been reduced by approximately 1000. 



